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I. INTRODUCTION 



The suggestion that the barrier distributions in subbarrier fusion reactions could be 
determined directly from the cross section data has led to a renewed experimental ac- 
tivity in the field [@-|^. Since barrier distributions are proportional to (P{Ea) / dE"^ , very 
accurate measurements of excitation functions at closely-spaced energies are required. Even 
with excellent data, smooth barrier distributions can only be obtained under certain model 
dependent assumptions [0 (i.e. - well-chosen energy spacing for calculating the second 
derivatives). Recently, it has been suggested that analyses using integrals over fusion data 
IQI provide model independent results, so these should be preferred over analyses that rely 
on differentiation of data. Specifically, the incomplete n-th moments of Ea, 

fJE) = n(n - 1) dE'iE - E'f'^E'aiE'), n > 2, (1.1) 

JO 

were proposed as an alternative for comparing model calculations with data. Unfortunately, 
these moments are not directly related to observables so it is difficult to give a physical 
interpretation for them. In addition, one of the main advantages of using barrier distributions 
is that they bring out important features in the data, but the the n-th moments in Eq. ( |1 . 1|) 
are even more featureless than the data from which they are calculated. 

The purpose of this paper is to point out how the average angular momentum, which 
is a measurable quantity §], is related to moments of the fusion cross section. In order to 
turn these relations into a practical tool to extract average angular momentum directly from 
fusion data, we study the angular momentum and deformation dependence of the effective 
barrier radius used in the formalism. In Section we review the relation between cross sec- 
tion and average angular momentum and introduce an expression for calculating the effective 
barrier radius. We derive improved relations in Section Q by using a better approxima- 
tion for the transmission probability for the £-th partial wave. This provides a satisfactory 
description of the relation between the cross section and average angular momentum for a 
spherical target. In Section |rV|, we present a qualitative discussion of deformation effects 
and propose a a simple way to include it in the effective radius. As a practical application 
of the method, we obtain the average angular momentum from the fusion cross section data 
in the ^^0+^^^Sm system and compare the results with the experiment. Finally, we draw 
conclusions from this work. 



II. THE GENERAL METHOD 

The idea of expressing average angular momenta in terms of integrals over functions 
of cross sections dates back to Ref. Due to absence of sufficiently good subbarrier 

fusion data, the full potential of this idea was not explored at that time. To introduce 
the concepts involved, we first review the general relations between moments of fusion cross 
sections and averages of powers of angular momenta in a one-dimensional barrier penetration 
picture. Although this is strictly valid only for spherical systems, it will provide the basis 
for extension to deformed systems. In the usual partial wave expansion, the total fusion 
cross section at the center-of-mass energy E is written as 
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£=0 

where the cross section for the i-th partial wave is 



(2.1) 



ME) 



7rr 
2^ 



{2i+l)TeiE) 



(2.2) 



T^{E) the transmission probability for that partial wave, and /i the reduced mass. For 
energies near the Coulomb barrier, one can approximate the ^-dependence in by using 



the s-wave penetrability at a shifted energy |Tl 



(2.3) 



where ixE?{E) is the effective moment of inertia of the system. The energy shift simply 
accounts for the change in the height of the barrier resulting from the centrifugal potential. 
Note that the effective radius is allowed to vary as a function of energy. In Section |T|, we 
will demonstrate how this approximation for the transmission probability can be improved. 



Substituting Eqs. ( p.2|) and (|2.3| ) into Eq. (|2.1|) , converting the sum over £ into an integral, 
and changing variables to 



E' = E- 



+ 1] 



2fiR^{E) 



(2.4) 



we obtain the expression 



a{E) 



7tR\E) rE 



E 



dE'To{E') 



(2.5) 



We will use Eq. (|2.5| ) to study the energy dependence of the effective radius from numerical 
calculations. 

The average angular momentum after fusion is assumed to be 



Following a procedure similar to the one used to obtain Eq. (|2 
momentum can be written as 



dE'ToiE') 



Ea{E) Jo 
Integrating by parts and using Eq. 



2fiR'^{E) 



2 {E-E') + - 



1/2 



(2.6) 

the average angular 



we arrive at the desired expression. 



h^Ea{E) Jo 



R^{E') 



^^{E-E) + - 



-1/2 



(2.7) 



(2.8) 
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which relates the average angular momentum to a moment of the fusion cross section. Note 



that the factor of 2 in Eq. (|2.8| ) is missing in Ref. [10|. For practical applications of Eq. ( p. 81) , 
it is important to note that, in order to obtain the detailed features of the average angular 
momentum reliably, the cross section should be interpolated rather than fit globally. We 
have found that a spline fit to the logarithm of the cross section works quite well for this 
purpose. This integral method with the assumption that the effective radius is constant. 



has been applied to some nearly spherical systems []r2|. The results are in agreement with 



the data within the experimental errors which, however, are rather large for a definitive 
confirmation that Eq. ( |2.8| ) with a constant radius works. The existing data for a variety 
of systems were also examined assuming a constant effective radius [|1^], but that analysis 



used a fit of the the cross section to the exponential of a polynomial so the features of the 
angular momenta are lost. 

Higher moments of the angular momentum can be found by following similar steps. For 
example, the second moment of the angular momentum is given by 

^^^^ + - ¥e^) L • ^2-^) 



A similar expression for {^'^) was given in Ref. [|1^, but {€) was neglected and R was assumed 



to be a constant. Under these assumptions, /2 in Eq. ( |Ll|) can be related to through 

{^')-^UE). (2.10) 

However, as has been pointed out previously [0,11^ , taking i? to be a constant is not a good 
approximation, especially for deformed nuclei. Therefore, Eq. ( p.lO| ) does not result in a 



reliable physical interpretation for /2. There is little experimental data for higher moments 
of i due to the difficulty of these measurements, so we do not pursue them here. The 
procedure for calculating them should be clear from the preceding discussion. 

In order to put Eq. ( |2.8| ) into use, we would like to understand the origin of the energy 
dependence in the effective radius better. For this purpose, we use Eq. ( p.5|) as the defining 
relation for R{E) and study its deviation from a constant value. We use values of a and 



To generated by the computer code IBMFUS [|T6[ which uses the interacting boson model 
(IBM) [|r^ to account for nuclear structure effects and evaluates transmission probabilities 
numerically in the WKB approximation. For fusion reactions involving rare-earth nuclei, 
IBMFUS has been shown to reproduce very well cross section and barrier distribution sys- 



tematics and average angular momentum data |T9|. It is important to note that the 
centrifugal energy is treated exactly in the WKB calculations (i.e. approximations such as 
Eq. (|2.3| ) are not used). In order to emphasize the effects of nuclear structure, we kept 
the masses of the projectile (^^O) and target (^^"^Sm) fixed and varied the quadrupole cou- 
pling parameter. Fig. |I] shows the results obtained for R{E) in three cases corresponding 
to the target nucleus being spherical (no coupling), vibrational (intermediate coupling) and 
deformed (strong coupling). These results demonstrate that the effective radius is not con- 
stant even for the spherical case and deviates more as the coupling to structure increases. 
The sharp increase in radius below the barrier with increasing deformation is obviously due 
to selective sampling of the longer nuclear axis. The origin of the energy dependence of 
radius in the spherical system is not that clear at this point, but the approximation used 
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in Eq. ( p.3|) is an obvious suspect. In view of the demonstrated energy dependence of the 
effective radius, extracting the average angular momentum from Eq. (|2.8| ) by assuming that 
R{E) is constant will not accurately predict {£) across a wide range of energies. 

Attempts have been made to parameterize R{E), but they have not been very successful. 
For example, a linear combination of the position of the barrier peak Rb and the Coulomb 
turning point 



Rc = ZiZ2e^/E . 



(2.11) 



has been suggested clS cl plausible choice O 



R{E) = 7]Rb + (1 - ii)R, 



c 



(2.12) 



In the spherical case, Eq. ( p.l2| ) provides a good description of R{E) in Fig. ^ with r] = 
0.78. However, is not clear how to include deformation effects in Eq. ( p.l2| ) in a physically 
meaningful way. Another expression for the effective radius, derived by assuming the nuclear 
potential to be an exponential tail in the region of the barrier, is given by [l^ 



R{E) = -Rc 



1 + (1 - 4a/i?, 



c) 



,1/2 



(2.13) 



where a is the nuclear surface diffuseness. This expression gives an energy dependence which 
is too strong for values of a in the range 0.6-1.2 fm. Also it doesn't make any allowance for 
inclusion of deformation effects, which are seen to have a significant influence on the shape 
of R{E). Clearly, a better understanding of R{E) is needed to make further progress. 



III. AN IMPROVED EXPRESSION FOR THE PENETRABILITY 

The prescription given in Eq. (|2.3D for approximating the i-wave penetrability by the 



s-wave penetrability at a shifted energy utilizes only the leading term in what is actually 
an infinite series expansion in A = i{i + 1). In this section, we derive the next term in this 
expansion and show the resulting corrections to the calculations presented in the Section 0. 
We also demonstrate how this can explain the part of the energy dependence of the effective 
radius which arises from the centrifugal potential. 

The effect of the angular momentum on the penetrability is usually taken into account 
by the shift it makes in the height of the potential barrier. The total potential for the i-th 
partial wave is given by 

Ve{r) = V^{r) + Vc{r) + J , ' , (3.1) 

where Vn and Vc are the nuclear and Coulomb potentials, respectively. Let denote the 
position of the peak of the i-wave barrier which satisfies 

dVf(r) , , 

= 0, (3.2) 



dr 

r=ri 



and 



5 



< 0, 



(3.3) 



r=ri 



then the height of the barrier is given by Vbi = Vi{ri). We make the ansatz that the barrier 
position can be written as an infinite series, 



rt = rQ + ciA + C2A^ H , 

where the q are constants. Expanding all functions in Eq. 
we find that the first coefficient is 



(3.4) 

consistently in powers of A, 



Cl 



where a is the curvature of the s-wave barrier 



a = — 



(3.5) 



(3.6) 



r=ro 



Substituting the leading order correction in the barrier position into Eq. (|3.1| ), we find 
that to second order in A the £-wave barrier height is given by 



Bl 



Vbo + 



4 a2 



rA 



(3.7) 



Therefore, an improved approximation for the ^-dependence in the penetrability is given by 



TeiE)=To\E- ^ 



(3.8) 



We give an alternative derivation of this expansion and discuss its validity in the Appendix. 

To examine the consequences of the improved expression for the penetrability, we repeat 
the steps outlined in Section P using Eq. (|3.8|) instead of Eq. ( p.3|) . To the leading order in 
1/a, this introduces the following correction to Eq. (|2] 



cr{E) = ^ ("dE'ToiE') 
E Jo 



arl 



[E - E') 



(3.9) 



Comparing Eq. (|2.5| ) with Eq. (|3.9|) , we find that the energy-dependent effective radius can 
be expressed as 



R\E) = rl 



4 !^dE'T,{E'){E-E') 



j^dE'T,{E') 



(3.10) 



This predicts the decrease in R{E) as the energy increases that was shown in Fig. 0. The 
calculation of the average angular momentum with the modified equations introduces a 
similar leading order correction. 



(0 



h'Ea{E) Jo 



dE' E'a{E') 



2^ 



{E -E') + - 



-1/2 
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{E - E') 



(3.11) 
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Since we have included a correction to re, this expression takes into account the shifts in 
both the height and position of the barrier for different partial waves. 

It is easy to show that the expression for effective radius given in Eq. ( |3.10| ) is consistent 
with the result obtained by weighted averaging over the barrier positions of all of the partial 
waves using Eq. (|3.4D, 

M = ?i^ = .„-^(A) + .... (3.12) 

To a first approximation, (A) is given by Eq. ( |2.9| ) with R{E) = tq. Substituting Eq. ( p.5| ) 
for E'a{E'), Eq. (|3l^ ) becomes 



in) = ro-^ r dE' r' dE"To{E"). (3.13) 



aEa Jo Jo 
After integration by parts, Eq. ( |3.13| ) becomes 

27rro 



(r,) =ro-^ f dE'To{E'){E - E'). (3.14) 
aEa Jo 

Substituting Eq. ( p.5|) with R{E) = vq for a and squaring the result, we recover an expression 
that is consistent with Eq. ( |3.1CI| ) to the first order in 1/a. 

The effect of the correction to the barrier position on the distribution of barriers is 



straightforward in the form given by Ackermann pO[, which is in terms of first derivatives 



of the angular momentum distribution. This expression can be written as 



ifi'R'E da,{E) 



where the shifted energy is 



Knowledge of the angular momentum distribution at an energy E allows the calculation 
of the barrier distribution between E — fi'^ tmax{^max + l)/2nR'^ and E, where imax is the 
largest angular momentum for which the partial cross section is measured. The first order 
correction to Eq. ( |3.15| ) is obtained by substituting the position of the £-wave barrier for the 
effective radius of the associated partial-wave cross section 

R^re^ro + ciA. (3.17) 

When the angular momentum distribution is determined at 65 MeV for the ^^0+^^^Sm 
system, the correction to D{E') is about 12% at 55 MeV and it will be less for higher 
energies. The other errors involved in calculating the barrier distribution are typically larger 
than that, so this correction can be safely ignored. 

In order to check whether or not the correction in Eq. ( |3.8D is sufficient to account for the 
shift in the peak of the barrier due to the angular momentum of the system, we calculated 
tq as a function of energy for a spherical target using Eq. (|3.9|). The same numerical values 
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of and Tq{E) obtained from IBMFUS as in Fig. |T] and the known value of a for the 
potential barrier were used in this calculation. The results for the s-wave barrier radius tq 
is shown in Fig. ^ and the effective radius R{E), extracted using Eq. ( p.5| ), is shown for 
comparison. The results for tq are nearly independent of energy as we expect. 

We would like to extract tq and a directly from cross section data. Eq. (|3^ ) can be 
rewritten as 

rE Aqr pE Ajr rE 

E(x{E) = 7crl dE'To{E') E dE'To{E') + — E'dE' To{E') . (3.18) 

JO a Jo a Jo 

Using partial integration, the two integrals can be expressed as 

„ dTo 



/ dE' UE') = EToiE) - - E'dE'^, (3.19) 
JO Jo dE' 



and 



^...'T,.0^^-rf..'^. (3.0, 

For E ^ Vbo, the right hand sides of Eqs. (p.l9|) and ( p.20| ) become E — qi and £'^/2 — g2, 
respectively, where 

l-E rIT 

q,= E'dE'^, (3.21) 



and 



dE' 



1 _,o dTo 



For energies above the barrier, dT^/dE' goes to 0, so qi and q2 become constants. Therefore, 
for high energies Ea{E) is a quadratic in E, 

Ea{E) = - (^) E' + (nrl + ^q,^E-- (nr^q^ + ^ ^2) • (3.23) 

Classically, 

"^^^ 6{E - Vbo) , (3.24) 



dE 



where Vbo is the barrier height, so qi = Vbo and q2 = V^q/2. Quantum mechanically, these 
expressions for qi and q2 are also approximately true for energies above the barrier. Using 
the values of Tq generated from the code IBMFUS [|T6| in Eqs. ( |3.21| ) and ( p.22| ), we have 
found that the error in these approximations at an energy of 70 MeV are both less than 
0.05% for an 0+Sm system with no coupling where Vbo ~ 59MeV. Therefore, the product 
of the cross section and the energy can be fit at high energies with the expression 

= nrl {E - Vbo) - — {E- Vbo)^ • (3.25) 
a 
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in order to determine a and tq. Of course, this requires high precision fusion data for energies 
above the barrier which may be difficult to obtain due to the competing processes. 

As a test of the formahsm, we apply the results derived in this section to the fusion cross 
section "data" generated by IBMFUS for a the system of ^^O + ^^^Sm, where both nuclei 
are taken to be spherical. First, a and tq are determined from Eq. ( |3.25| ) by a fit to the a 
"data" . Then these values are employed in Eq. and the average angular momenta are 

extracted from the a "data" through numerical integration. Fig. ^ shows a comparison of 
the average angular momenta calculated using Eq. with those obtained from IBMFUS 

directly. The agreement is very good at all energies. For reference, results obtained from 
Eq. ( p.8| ) assuming a constant radius are also shown (dashed line) . As expected, the modified 
expression ( |3.11|) leads to a clear improvement at high energies where the effective radius 
varies most (cf. Fig. |^). This example gives confidence that one can use Eq. ( p.ll| ) in 
extracting average angular momenta directly from the fusion data for spherical systems. 



IV. TARGET DEFORMATION EFFECTS 



Now that we have an improved description of the effects due to the centrifugal barriers, 
we consider the effects of target deformation. The shape of an axially symmetric deformed 
nucleus can be described by 



7^t(^^) = 7^io(l + /3r2o(cose)) 



(4.1) 



where TZto is the radius for an undeformed target. In order to qualitatively study deforma- 
tion effects, we approximate the target with a simple two-level system which displays the 
important features. As shown in Ref. [^, fusion of a deformed nucleus with finite number 
of levels [n) can be described by sampling n orientations of Eq. ([4.1|) with their respective 
weights. For a two-level system, the orientations 6i = 70.12° and 62 = 30.55° contribute 
with the weight factors wi = 0.652 and W2 = 0.348, respectively. 

To proceed, we need to find out how the barrier position and height changes with orien- 
tation, which can be calculated from the total potential in a straightforward manner. For 
the nuclear potential, we use the usual Woods-Saxon form with the target radius given by 
Eq. (O) 



VN{r,9) = -Vr 



NO 



1 + exp 



-1 



(4.2) 



The Coulomb potential can be calculated from a multipole expansion and, to leading order 
in P, is given by 



Voir, 



ZiZ2e' , 3/37^2oy2o(cos( 



A / ^ 3^to(^i(^l^^^ 



(4.3) 



where A = Z\Z2e^ ■ Using the notation of Section p| , the equation for finding the peak of 
the s-wave barrier is 



a [l + exp((r-7^p-7^^)/a)]^ 



0, 



(4.4) 



r=ro 
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and the height of the s-wave barrier is 



A 



to 



V, 



NO 



1 + exp 



ro-TZp- Tit 



(4.5) 



The 6'-dependence is suppressed in the above equations for convenience, but both ro and 
Vbo depend on the target orientation. 

The rate of change in the barrier height due to the deformation is given by 



dVBo dVso dro _^ ^ATZto 



dlZt 



5r^ 



Vm exp ((rp - Tip- Tit) / a) 
a [1 + exp ((ro -Up- Tit) /a)] 



2 • 



(4.6) 



In Eq. (|4.6| ), dVBo/dro = by definition and using Eq. (|4^ ), the second term can be 
simphfied to give 



dVr 



BO 



dTZt 



nt=nto 



To find a similar expression for the barrier position, we differentiate Eq. 
to Tit 



(4.7) 
with respect 



2 A 36Anto{nt~nto) 

■3 ~r - c; 



dro QATZto ( dro 



dn 



5r^ 



X 



Vno exp ((rp - Tip- Tit) / a) 
[l+exp((ro-7^p-7^^)/a)]' 



\dnt J 

2 exp {{ro-TZp -Tit) /a) 
[l + exp((ro-7^p-7^^)/a)] 



0. 



(4.8) 



Using Eq. ( [4.4| ), we can solve for the rate of change in the s-wave barrier position due to the 
change in the target radius 



drn 



where 



d-JZt 



Q 



7^t=7^^o 



Q - 9a7^ffl/5rg 
Q + 2a/ro 



^ (--Vso]-l. 



(4.9) 



(4.10) 



Va^o Vrp 

For a given orientation 6, the shift in the s-wave barrier height is given approximately by 



6V, 



dV, 



BO 



BO 



dn. 



A SATZ 



to 




pntoPi{cost 



(4.11) 



and, similarly, the shift in the s-wave barrier position is approximately 



Srn 



dro 



dn. 



dro 



nt=Tito 



dn, 



7^^=7^to 




/37^^pP^(cos( 



(4.12) 



These expressions account for the changes due to deformation fairly accurately as can be 
seen in Fig. H. 
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The total fusion cross section for a two-level system is given by 



ariE) = w^a{E, Ai) + W2a{E, A2) , (4.13) 

where (j{E, Aj) is the cross section for the i-th orientation and Aj = P2{cos 9i). To simplify 
the notation, we introduce 



5 dro 



7^^o, (4.14) 

nt=nto 



47r dTZt 

so that, for a given orientation, the peak of the s-wave barrier in Eq. ( |3.9| ) is replaced by 

ro^ro + F/3X,. (4.15) 
After this substitution, the cross section for each level becomes 

a{E,X,) = I J^'^dE' {rln{E\\i) + F(3r^Kn{E\K)+F^(3^\lT^{E\\i)] 

_il f''dE'To{E',X,){E-E'), (4.16) 
ah Jo 

where To{E,Xi) is the transmission probability for i-th orientation. The curvature of the 
barrier, a, is expected to have a second order dependence on hence it is assumed to 
be constant in this leading order calculation. Defining the coupled s-wave transmission 
probability as 

To^(E) = wiToiE, Ai) + W2To{E, A2) , (4.17) 
the total cross section becomes 



/(f dE' [wAME', Ai) + W2X2To{E', A2)] 



aT{E) = ^{Tl + 2FpT, 



, /o^ dE' [w^\lToiE', AQ + W2\lToiE', A^)] ] -^.^.c.^n 
S^dE'T^{E') jJo ' 

dE'T^{E'){E-E'). (4.18) 

Jo 



~aE 

The cross section can be written in the form of Eq. (|2.5| ) as 

aAE) = ^RUE) j^E'T^iE'), (4.19) 
so that the effective radius with coupling is 

RUE) = rl + 2F,}r, Jo^ N.A.r„(g. AQ + ..A,r„(g. A,)] 

' ° Jo^dE'T^{E') 

^p2^2 lodE' [wi\lToiE',\i)+W2\lTo{E',X2)] 



J.'^dE'T.^iE') 

4 J,'=dE'T,^{E'){E-E') 
a SodE'T§{E') 
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(4.20) 



Since wiXi = W2X2 and both To{E', Ai) and Tq{E' , A2) approach one very quickly for energies 
above the barriers, the second term in Eq. ( 4.201 ) becomes zero for high energies. On the other 
hand, the third term of that equation is always positive, so by comparison with Eq. ( p.lO| ) 
it is easy to see that the effective radius is slightly higher at large energies for the deformed 
case than for the uncoupled case. This argument also holds for multi-level systems, since 



j Pi{cos9) d{cos9) = 



(4.21) 



The effective radius predicted by Eq. ( [4.20|) (dashed curve) and the result from definition 
of Eq. ( |4.19| ) (solid curve) are shown in Fig. ^. That the two curves are in good agreement 
is an indication that the curvature is approximately unchanged as assumed in Eq. ( [4.201) . 
This simple model exhibits the main features seen in Fig. 0; at low energies the difference 
between the deformed and spherical cases becomes larger and this effect increases with the 
deformation. 

The preceding argument provides an explanation of how the effective radius varies with 
the energy. Unfortunately, it is not easy to incorporate the deformation effects in the 
calculation of average angular momentum as was done for the centrifugal barrier in Section 
lT|. To make progress, we take a phenomenological approach and introduce the quantity 



pUe) 



Eg/ 71 + (4/a) /^f dE' To{E'){E - E') 
lo^dE'ToiE') 



(4.22) 



For a single barrier, po{E) is simply the location of the s-wave barrier ro (cf. Eq. ( |3.9|) ), so it 
is actually energy independent. When couplings to target deformation are introduced, there 
is a distribution of barriers. In this case, pl{E) is a suitable average of the location of the 
barrier peaks. Using the values of cr and To generated by IBMFUS in Eq. ( [4. 22]) again, we 
calculate Po{E) for the three quadrupole coupling strengths used in the previous section. In 
contrast to the effective radii in Fig. |l|, the results shown in Fig. ^ level off at high energies. 
They can be parametrized using a simple Fermi function 

1 + exp((E - Vbo)/W) 

Here tq is the asymptotic value at large E and Vgo is the barrier height, which are determined 
from the fusion data using Eq. ( ^.25| ). + 5 corresponds to the asymptotic value at low 
energies and can be calculated from Eq. ( |4.1| ) at = 0. The only quantity in Eq. (f4.23|) 
that is not determined from data is the width W . The fits to the curves in Fig. ^ results 
in values around 2 ± 0.3 for W with a mild dependence on (3 (increases with 13). Since the 
precise value of W makes no tangible difference, the slight uncertainty in its value is not 
important for the purposes of this paper. 

The averaged radius of Eq. |4.23| can now be used in Eq. (|3.11|) in the place of ro to obtain 
an improved expression for average angular momentum. 



h^Ea{E) 



'dE'^ 

pI{E') 



-1/2 



1 - 
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apliE) 



(E - E') 




a \pI{E') pI{E)^ 



-^^{E-E) + - 



(4.24) 
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Note that when po{E) = tq, the last term vanishes and the above expression reduces 
to Eq. ( p. 11 ). The average angular momentum extracted from the fusion "data" using 



Eq. ( [4.24| ) is compared to the {£) values obtained from the same IBMFUS calculation in 
Fig. 1^. The agreement is very good at all energies. In contrast, the constant radius results 
underestimate {£) by about one h. 

To demonstrate the utility of this method in extracting {£), we apply it to the ^^0+^^^Sm 
system for which quality cross section data exist 0. In Fig. we compare the {£) values 
obtained from Eq. ( |4.24|) and IBMFUS with the experimental data |2^. The two calculations 
are consistent with each other but slightly underpredict the experimental values. 



V. CONCLUSIONS 



We have described two major reasons for the energy dependence of the effective radius. 
The effects of the centrifugal barriers are described by using an improved approximation for 
the penetration probability. This also leads to a better relation between the cross section 
and the average angular momentum. The effects of target deformation are described with 
a simple model which reproduces the features of the effective radius. Finally, we present a 
phenomenological expression for the position of the s-wave barrier as a function of energy 
and show how it can be used in extracting the average angular momentum from the fusion 
cross section. Comparison of this method with numerical calculations shows that it predicts 
the average angular momentum from the fusion data reliably, and hence it can be used as a 
consistency check in cases where quality data are available for both quantities. 
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APPENDIX: THE VALIDITY OF THE EXPANSION OF THE PENETRABILITY 



In this appendix, we discuss the validity of Eq. ( |3.8| ) which approximates the transmission 
probability as a power series in A = £{£ + 1). We will do this by using the linearized forms 
of the WKB penetration integrals. The penetration probability of an ^-wave through a 
one-dimensional barrier is given by 

T,{E) = [l+exp{2Se{E))]-\ (Al) 



where the WKB penetration integral is 



V -^rii 



Voir) + 



h^£{£ + 1) 
2/ir2 



E 



1/2 



(A2) 



Using Abelian integrals, one can show that for energies below the barrier 
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TT 2fi 

VE'-E " 2 V F 



r2l 



dr 



Voir) + 



+ 1] 



-E 



(A3) 



where Vbi is the height of the £-wave potential, Vo(r) is the s-wave barrier, and ru, are 
the turning points of the £-wave barrier for energy E. In Refs. [pT| , p3| , the energy derivative 
of Eq. (^, 



dS,{E')/dE' 
? VE'~E 



TT 



!2fi 



was used to find the barrier thickness. One can also take the derivative of Eq. 
respect to A to obtain another useful identity 



dSe{E')/dA 



TT 



1 

r2e 



1 



(A4) 
with 

(A5) 



y/E' -E 

The last two equations can be used to check the consistency of an expansion of the form 

S,{E) = So{E - eiA - esA^ - ...) , (A6) 
are independent of the energy. Substituting the relationship between the 



where 61,62 
derivatives, 



dSe{E) 



(ei + 2e2A + ... 



dS.iE) 



into Eqs. 



and 



OA ^ ' ^ ' dE 

we find that the assumption in Eq. 

(ei + 2e2A + ...) 



is consistent if 



'2fxrur2e 



(A7) 



(A8) 



The right hand side of this equation has a slight dependence on energy (through the energy- 
dependence of the turning points) whereas the left hand side is independent of energy in our 
approximation. Our analysis in Section |lll| is equivalent to approximating the right hand 
side of Eq. (|7^) by a constant, i.e. 



rur2i ~ r^ , 



(A9) 



where is the position of the peak of the £-wave barrier. This is a very good approximation 
at energies near the barrier height, but the error increases as the energy gets lower. However, 
for the ^^O +^^^Sm system, even at 7 MeV below the barrier (lower than the fusion cross 
section has yet been measured) the error in this approximation is less than 4 % for the bare 
potential used in Ref. [Q. Upon setting £ = in Eq. (|A9|) , we recover 



ei 



2 ' 



(AlO) 



in agreement with Eq. 



The higher-order terms in A also agree with the calculations 



in Section |T|, so we are justified in expressing Ti{E) as a power series in A. 



14 



It should be noted that although A is not a small parameter, there is a natural cutoff 
Acr in this parameter given by the condition 



For values of A greater than this, fusion will not occur. Due to this cutoff, the second term 
in Eq. (|3.8| ) will always be larger than the third term. Therefore, the additional term can 
be considered a correction. 
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FIG. 1. The effective radius R{E) extracted from fusion calculations for the sys- 
tem using Eq. (p.5|). The curves correspond to spherical, vibrational and deformed nuclei with 
quadrupole coupling strengths V2 = 0, 0.13 and 0.26, respectively. 
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FIG. 2. The s-wave barrier position tq (solid curve) extracted from fusion calculations for the 
system with no coupling using Eq. ( |3.9D . For comparison, the values of R{E) (dashed 



curve) calculated for the same reaction using Eq. (2^) are also shown 
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FIG. 3. Comparison of average angular momenta calculated using Eq. ( 3.11| ) (solid curve) and 
Eq. (|2.8| ) (dashed curve) to the values calculated by IBMFUS (points) for the ^^0+^^^Sm system 
with no coupling. 
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FIG. 4. The solid curve is the total potential for the system when the target is taken 

to be spherical. The dashed (A = —0.327) and dot-dashed (A = 0.613) curves are the potentials 
for two-level approximation with (3 = 0.25. The arrows show the shifts predicted for the barrier 
peaks by Eqs. pH] ) and ( p^ . 
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FIG. 5. The effective radius from a two-level calculation for the system with 

P = 0.25. The solid curve is calculated using the definition, Eq. ( [4.19|) , and the dot-dashed curve 
usmff Eq. (|]2|), which assumes a is constant. For comparison, the dashed curve is the result for 
a spherical target. 
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FIG. 7. Comparison of average angular momenta calculated using Eq. ( |3.11 ) with a constant 
value (dashed curve) and Eq. ( p3| ) (sohd curve) for po to the values calculated by IBMFUS (points) 
for the ^^0+^^^Sm system with the quadrupole coupling strength V2 = 0.26. The parameters used 
were ro = 10.6 fm, 5 = 1.22 fm, Vbo = 59 MeV, and W = 2.3 MeV. 
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FIG. 8. Average angular momenta data y] for the system 0+ Sm. The soHd curve shows 



the results using Eq. (4.24) and dashed curve those from IBMFUS ||T 
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